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Abstract 

We extend a theory of first order h corrections to Gutzwiller's trace formula for 
systems with a smooth potential to systems with discrete symmetries and, as an 
example, apply the method to the two-dimensional hydrogen atom in a uniform 
magnetic field. We exploit the C4t,-symmetry of the system in the calculation of the 
correction terms. The numerical results for the semiclassical values will be compared 
with values extracted from exact quantum mechanical calculations. The compari- 
son shows an excellent agreement and demonstrates the power of the h expansion 
method. 
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1 Introduction 



Semiclassical theories have become very important for a deeper understanding 
of quantum systems, and Gutzwiller's trace formula [1] has become a powerful 
tool for classically chaotic systems. It provides a semiclassical approximation of 
the quantum level density in terms of classical periodic orbits. In a systematic 
expansion of the level density in powers of h it can be considered as the leading 
order. Higher orders of this asymptotic expansion have been developed in 
several studies [2,3,4], but for a long time were only tested for billiard systems, 
i.e., systems with hard walls instead of smooth potentials. By extending an 
expansion which was derived by Gaspard et al. [3,4], Gremaud [5] developed h 
corrections to Gutzwiller's trace formula for quantum systems with a smooth 
potential. He presented a first-order h correction term to Gutzwiller's trace 
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formula and obtained numerical results for the diamagnetic hydrogen atom. 
Comparisons with values extracted from exact quantum calculations showed 
that the semiclassical results can be computed with very high accuracy. 

However, in the theory presented in [5] some important topics have not yet 
been considered. For example, the H corrections cannot be calculated for clas- 
sical orbits which have a turning point (i.e., the velocity vanishes completely 
at this point). Furthermore, discrete and continuous symmetries have not been 
included. 

If there are discrete symmetries the eigenstates of the quantum system split 
up into several subspaces. In these subspaces, classical orbits which are not 
periodic without a symmetry transformation contribute to the level density. 
Preliminary results for the diamagnetic hydrogen atom, which is an exam- 
ple system with discrete symmetries, were published in [6] without explain- 
ing which modifications in the numerical calculations the introduction of dis- 
crete symmetries entails. In particular, an analysis of the differential equations 
which have to be solved to obtain the correction terms is not presented. In this 
paper we want to give a transparent derivation of the modifications necessary 
for the calculation of the correction terms in systems with discrete symmetries. 
We will show in detail how such symmetries of the Hamiltonian have to be 
taken into account. It will be explained in which equations it is necessary to 
introduce symmetry operations and which transformations lead to the correct 
boundary conditions of the classical Greens function, which is an essential part 
of the correction terms. Furthermore, we want to give a deeper understanding 
of the significance of the symmetry of the system for the quantum spectrum 
and the connection to periodic orbits. 

We will apply the method to the diamagnetic hydrogen atom. The three- 
dimensional hydrogen atom in a uniform magnetic field has a continuous 
symmetry, namely the rotational invariance around the magnetic field axis. 
Continuous symmetries have a substantial influence on the correction terms. 
Actually one cannot obtain the correct results for the first-order h correction 
for the three-dimensional hydrogen atom with the formulas presented because 
the rotational invariance leads to additional contributions to the first-order 
corrections. Therefore, we will consider the hydrogen atom as a pure two- 
dimensional system as it was done before in [5]. If one uses semiparabolic 
coordinates, the potential of the two-dimensional diamagnetic hydrogen atom 
exhibits a discrete C4^-symmetry. The h corrections will be calculated for 
selected periodic orbits taking into account that discrete symmetry. The semi- 
classical results will be compared with the analysis of exact quantum calcu- 
lations. The agreement between the results of both methods turns out to be 
very good. 

The outline of the paper is as follows. In section 2 we will first give a summary 
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of the derivation of the H corrections without taking into account discrete 
symmetries. In section 3, wc will introduce the hydrogen atom in a uniform 
magnetic field with all aspects relevant for the calculation of the correction 
terms. We will also discuss the influence of the rotational invariance of the 
three-dimensional hydrogen atom. Then we will extend the h corrections to 
discrete symmetries, calculate results for the two-dimensional hydrogen atom 
and compare them with exact quantum mechanical calculations in section 4. 



2 Semiclassical approximation of the quantum level density 



As mentioned in the introduction the h corrections to Gutzwiller's trace for- 
mula for systems with a smooth potential were derived by Gremaud [5] based 
on earlier work by Gaspard et al. [3,4]. In this section, for the reader's conve- 
nience we will briefly summarize all major steps of this derivation by following 
the way described in [5] in order to be able to include discrete symmetries in 
section 4. 



2. 1 The trace of the propagator and the first part of the H correction 



The starting point for the derivation of the semiclassical level density which 
is used in this theory is a discrete version of Feynman's path integral repre- 
sentation of the propagator 



K{q, qo, T) = VTW / ^gi dgs ■ ■ ■ dg7v_i 



{2mhAt)- 

X exp 



for a system with / degrees of freedom and a time independent Hamiltonian: 

H = ^ + V{q). (2) 

K{q, qo, T) represents the propagation of a particle going from % to q^ = q 
in time T with At — T/N. L{q, q) denotes the classical Lagrangian of the 
system. 

For the calculation of the level density one has to know the trace of the prop- 
agator, which can be obtained by setting % — g in (1) and by integrating over 
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this variable: 



K{T) = [ K{qo, qo, T) dg'o = ^ — w I ^qo dqi d?2 • • • d^jv-i 

^ (3) 



(27ri;iAt) 

X exp 



After expanding the exponential function and including all terms which pro- 
duce contributions to the leading order in h and to the first-order correction, 
the integral is evaluated as described in [4] . The result in the hmit At — > is 
given by: 



K{T) = Y: Kf\T) |l + ^ j^' dto Cu(T, to) + O {h') \ . (4) 



In this formula the sum runs over all periodic orbits of the classical system 
described by the Hamiltonian (2) with period T. Tp represents the time period 
of the primitive periodic orbit, i.e., if the orbit consists of multiple repetitions 
of a periodic orbit, one considers only the basic traversal in this case. Kf\T) 
stands for the leading order of the semiclassical approximation. It is identical 
with Gutzwiller's [1] result, and is given by: 



Kf\T) - ^ ^ 



V27fn J\dETdei (me(T) - 1) 



X exp 



(5) 



In this expression one can find the (2/ — 2) x (2/ — 2)-dimcnsional stability 
matrix m, which represents the stability properties of the periodic orbit. The 
determinant det {me{T) — 1) is related to the 2/ x 2/-dimensional monodromy 
matrix. This relation is shown in section 2.4.2. A further term is the classical 
action 

wt'\T)= r L{qA)dt, (6) 
Jo 

/li is the Maslov index, and OeT is the derivative of the time period with 
respect to the energy E. In equation (4) Ci e is the abbreviation of the first- 
order correction in h. As one can see in (4) one has to use the average of 
Ci^(T, to), where to parametrizes the periodic orbit. The integral over to runs 
from to Tp, i.e., one has to integrate over the whole primitive periodic orbit. 
The explicit expression of Ci^(T, to) needs lengthy calculations, which are 
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discussed in [3,4,5] in detail. The result reads: 

1 /"-^ 

Cie{T, to) ^ g Kni2j3M (^)) Sin2{'t, t'l to)Si3i4{t, t; to) 



24 Jo 

^ ^i3ji(^> 'to)Gj2j3i'^' ^ ^''j ^o) + 2^jjjj(t, t'; to)Qi2j2{t, t'; to)^i3j3(^> ^'j 



dil<M2i3 ^in(0,i;io)^i.i3(^,^;^o), (7) 



where 

are derivatives of the potential evaluated at the point g^*^^^ (i) on the classical 
orbit. Qij are the components of the classical Green's function, which is a 
solution of the linearized equation of motion: 

(-4 "Sf 5<'-'''^i^<'-'')- 

Due to the factor l/jg^'^'^p in equation (7), the correction term is singular for 
vanishing velocities, and because of the integral over to in equation (4), the H 
correction term diverges for orbits with a turning point. 



2.2 The classical Green's function for the trace of the propagator 



The classical Green's function is an essential part of the first-order h correction 

to the trace of the propagator. In order to be able to look at the symmetry 
transformations of Cie in equation (7), one has to know its structure. In this 
section we discuss all important parts of this structure. 

Prom the derivation of the correction term Cn, it follows that the classical 
Green's function has to fulfil the boundary conditions [5] 



g{o,t') = giT,t'), 
VtMo,t')^VtMT,t')^o, 



(10) 



with the projection operator Vto along the direction of the classical orbit at 
time which has the form 



q {to) (g) q{t o) 
q{to) ' 



qMqjito) 
q{to) ' 
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furthermore, Qt^ = 1 — Vt^, with the /-dimensional unity matrix 1. If we use 
the notation 

g_(t, t') = g{t, t') for < t < t', 

g^{t,t') = g{t,t') for t'<t<T, 

it is possible to write the classical Green's function as a product [5] 



(12) 



^M{t) 



B±{t') 



(13) 



with the 2/ x 2/ monodromy matrix M{t) and the four / x / matrices A± and 
B±. The monodromy matrix is a symplectic matrix, which can be obtained 
by solving the linearized Hamiltonian equations of motion 



M{t,T) 



'dXdX 



^M{t,T). 



(14) 



Here, S is the matrix 



/ ol^ 



V 



1 



and 1 is the / x / unity matrix. For a Hamiltonian of the form (2) it has the 
structure: 



/ 



M{t) 



\ 



(15) 



J2(t) Ji{t) 

where Ji and J2 are f x f matrices. The boundary conditions at time t = t' 



are 




(16) 



Exploiting the formulation (13) of the classical Green's function and the sym- 
plecticity of the monodromy matrix M{t), one can formulate the boundary 
condition (16) as follows: 



A4t') 



A.{t') 



\B_{t') 



-M{t 




(17) 



If this expression is combined with the condition 



1 

\o QtoJ 



A_{t') 
B-{t') 




M{T) 




(18) 
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which is identical with the first and third condition of (10), it leads to the 
matrix equation: 




[M{T) 




M{T) 



(19) 



This equation has to be solved in order to obtain the matrices A_ and B_. 
When A_ and B_ are known one has the solution for the classical Green's 
function (13). Unfortunately, the matrix 




[M{T) 



is singular but, as was shown in [5], there is a solution for classical trajectories 
which have no turning point (i.e. the velocity q{t) never vanishes completely 
along the orbit), viz. 



— X — Xq + 



q{to) ® g(to) 



(20) 



where Xq represents a particular solution of equation (19), which can be deter- 
mined by a singular value decomposition (see e.g. [7]). The second condition 
of equation (10) has been used to derive this result. 

In summary, it is possible to obtain the classical Green's function by solving 
the linearized equation of motion (14) of the monodromy matrix and by cal- 
culating A_ and B_ from equation (20). Then the correction term Cie{T,to) 
follows from equation (7). In equation (20), the factor l/|g('^')p appears in 
the construction of Ci^(T, to) a second time. The reason for this lies in the 
application of the projection operator Vto in the second condition of (10). 

The integrals and double integrals in equation (7) can be transformed into a 
set of ordinary differential equations, which can be computed effectively. The 
transformation of the integrals is discussed in [5] in detail. 



2.3 The trace of the quantum Green's function and the second part of the H 
corrections 



The next step on the way to the level density is the trace of the quantum 
Green's function, which can be obtained from the trace of the propagator via 
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the semi-sided Fourier transform 



1 f°° 

G{E) = J2Ge{E) = j:- dTexp 

/> 1''' "^0 



-ET 



h 



Kt{T). 



(21) 



One has to include all contributions to the leading order as well as the first- 
order h correction in the semiclassical approximation of the integral. The meth- 
ods which are used to evaluate this integral are discussed in [4]. We do not 
want to repeat the calculation but present the results obtained in [4] and [5], 
which can be summarized as follows: 



G^{E) 



1 



: exp 



-S, (To)--/x, 



i^^|det(m,(ro)-l), 

x{l + ih (CieiTo) + Cfr^(ro)) + O (/l^) } , (22) 



where the leading order, belonging to the 1 in the curly brackets, is known 
from Gutzwiller's trace formula [1] . The reduced action 



Sf\To)^Wr{To) + ETo 



(cl). 



(23) 



is used in this expression, and (le represents the Maslov index. The first-order 
h correction consists of two terms. The first. 



1 /■'^op 

Cu(To) = — Cie{To,to) dto, 

I Op Jo 



(24) 



is the same as in equation (7) but now with fixed period Tq, and the second 
is the contribution from integral (21) to the h correction, and is given by: 



2WP{To 



^(2) /rp 



K3), 



+ 



2 W^iTo)^ ^WPiTo)^ 2iwP{Tof 



(25) 



and are the derivatives 



.in) 



wt\T) 



d"]]f^iT) 



Qrpn 



of the action W^'^^^ (T) and of the logarithm of the amphtude of the trace of 
the propagator 



Coi{T) = In 



T 



OeT det {rntiT) 



(26) 



respectively. 
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2-4 Contributions to the correction terms of the trace of the semiclassical 
Green's function 



In this section we want to summarize the contributions to the second correction 
term Cj^^ which are important for the understanding of the symmetry prop- 
erties of this term. As we will see, all parts of equation (25) can be obtained 
as numerical solutions of differential equations of some "new" coordinates. 



24.1 Derivatives ofWf\T) 
Because of the relation 



dT 



-E{To)^-H{X{t,To)) 



(27) 



T=To 



all higher derivatives of the action W^'^^^ (Tq) with respect to the time period of 
the orbit can be expressed as derivatives of the classical Hamiltonian, where 
X{t,To) = (g(t, To),p(t, To)) represents the phase space vector of the orbit. 
Introducing the derivatives 



^ ^*'^o^ - — — 



(28) 



T=To 



which appear as coefficients of the Taylor expansion 

00 1 

X{t,T) = X{t,To + ST)^Y1 ^X^''\t,To){STy 



n=0 



nl 



(29) 



of the orbit X{t, T) around a reference orbit X{t, Tq) with period Tq, the higher 
derivatives of w/''^(ro) are: 

wPin) = -H,,{x{t,n))x^'\t,n), 

wPin) = -H,ij{x{t,n))xi'\t,n)x^'\t,n) - H,,{x{t,n))x\^\t,T,), 



W^f (To) = -H,,k{X{t,T,))Xr\t,n)Xf\t,n)Xl'\t,T,) 



(1)/ 



'(1)/ 



'(1), 



- m,i^{X{t,%))xf'\t,%)Xf\t,%) - Hi{X{t,To))X^'>{t,To) 



'(2) 



(3)/ 



(30) 



From Hamilton's equations of motion for the phase space coordinates 



x(i,ro) = s 



dH{X{t,To)) 
dX 
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one can infer the equations which govern the motion of the derivatives X^") (i, Tq) : 
X^'\t,To) = J:,jHjk{X{t,To))xi'\t,To), 

(i,To) = ^ijHjki{x{t,n))xi'\t,n)xi'\t,n) 

+ j:,jHj,iXit,To))xi'\t,To), 
X^'\t, To) = J:,jHjMm{^it, To))xi'\t, To)X^'\t, T,)X^^\t, T,) 
+ 3E,,//,,-fe,(X(t, n))X^^\t, To)xj:^\t, To) 

+ E,,i/,,fc(X(t,To))xf (t,To). 

These equations are inhomogeneous differential equations. Their solutions are 
of the type 

x«(i,ro) = M(t,ro)^«(o,ro), 

X^'\t,To) = M(t,ro)X(^)(0,To) + F(^)(t,To), (32) 
X^^\t,To) = M(t,To)X(3)(0,ro) + F(3)(t,To), 

where iVf(i, To) is the monodromy matrix, which was introduced in section 
2.2. Together with the condition 

X{t^O,To)^X{t = To,To) (33) 

for periodic orbits, the solutions (32) lead to the following equations for the 
initial values: 



[1 - M{To, To)] X« (0, To) = X{To, To), 

(34) 



[1 -M(To,To)]X(2)(0,To) = 2X«(To,To) +X(To,To) + F(2)(To,To) 



[1 -M(To,To)]X(3)(0,To) = 3X(2)(To,To) +3X«(To,To) +X(To,To) 

+ F(3)(To,To). 

The matrix [1 — M{To,To)] is singular. A particular solution can be found 
via a singular value decomposition. For a general unstable orbit the kernel is 
one-dimensional, i.e., the solution space is in general one-dimensional (see [5]) 

(2) 

but, as also shown in [5], IV/ (To) can directly be computed if a particular 
solution Xq^ (0, To) from (34) is known: 

wPiTo) = -VH{X{To,To))-X^'\q,To). (35) 

Similar relations can be found for W^/^'*(To) and W^'^\To), in which also only 
a particular solution of equation (34) is required. 
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Derivatives ofCoe{T) 



If one looks at the definition ofCoe{T) in equation (26), one can easily convince 
oneself that 



l_ l wf\To) 1 d 
To^2wP{T,) 2 dT 

K4) 



ln|det (m(T) - 1)| 



T=To 



1 lWP{To) 1 / W^f^Tc 



^2 



1 d^ 



(36) 



2 dr2 



ln|det (m(T) - 1)| 



T=To 



The derivatives of the action are already known from section 2.4.1. In the next 
step one has to find an expression for the derivatives of In |det {m{T) — 1)|. 
Following again the method presented in [5] wc can introduce a new matrix 
N{T) whose determinant is identical with det {rn{T) — 1). It is given by 



where 



N{T) = M{T) - (l - V\\{T) - V±{T)) , 
X{0,T)^X{0,T) 



P||(T) = en ® eii 



X{0,T) 



(37) 



(3J 



and 



X{0,T) 

are the projection operators parallel to the flow of the classical orbit and 
perpendicular to the energy shell, respectively. These relations lead to the 
expressions 



— ln|det(m(r)-l) 



Tr N{To) 



-1 diV(T) 



T=To 



dT 



(40) 



T=To/ 



dT2 



ln|det {me{T) - 1) 



N{To) 



= TriN{To) 

T=To V 

-1 dN{T) 



.1 d''N{T) 



dT 



dr2 

N{To) 



T=To 



rp rp 

-1 dN{T) 
dT 



(41) 



T=Toy 



with 
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dN{T) 

dT 

dT2 



dM{T,T) ^ dV\\{T) 



dT 



dT 



d^M{T,T) , d2P||(T) 



+ 



dT2 



dT 
dT2 



(42) 
(43) 



By using the derivatives X^^\t,To) and the definition of the projection oper- 
ator Vn (T) , one obtains: 



dViijT) 
dT 



d-V\\{T) 

dr2 



X(^)(0,T) ®X(0,T) +X(0,T) ®XW(0,T) 



X{0,T) 

i(o,r) •x(i)(o,r) 



- 2 



X(0,T) 



^11 (7^), 
X{0,T) 

-2 fx(2)(0,T) ®X(0,T) +X(0,T) ®X(2)(0,T) 



+ 2xW(o,r)(g)XW(o,r)) + 8 



x(o,r) •xW(o,r) 



- 2 



X(0,T) 

X(0,T) •X(2)(o^2^) X(i)(0,T) •XW(0,T) 



X(0,T) 



-2- 



X(0,T) 



x(o,T)-xW(o,r).^. 



x(o,r) 

+ i(0,T)®i«(0,T)). 



(X«(0,T) ®X(0,T) 



The last pieces which are needed for calculating Coi{T) are the first and 
the second derivative of the monodromy matrix M[T,T) with respect to the 
period T: 



dM(T, T) _ dM{t = T, T) dt dM(t = T, T) 
dT ~ di dT ^ ffT 

= M(T, T) + M^^\T,T), 
d^M{T T) ^ ^ ^ ^^(1) ^ ^(2) ^ 

dT2 



(44) 



Similar to the procedure which was used in section 2.4.1 to arrive at the 
equations of motion for the X^"'\t, Tq), it is possible to obtain the differential 
equations governing the evolution of the derivatives M^") (T, T) from differen- 
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tial equation (14) of the monodromy matrix. The equations read 

M^}^(t,T) = {H,,U^it,T))Mijit,T)X^\t,T) + H,i{X{t,T))Mil\t,T)) , 

+ 2H^kUX{t, T))Mif it, T)X« (t, T) 

+ H,um{X{t, T))Mij{t, T)X^^^ {t, T) + H,ki{X{t, T))Mlf {t, T)) , 

(45) 

with initial values 

M(")(i = 0,T) = 0. (46) 

Now we have obtained all results necessary for the calculation of the first- 
order h corrections to Gutzwiller's trace formula without taking into account 
symmetries. In the next sections we will apply this theory to the diamagnetic 
hydrogen atom. 



3 The hydrogen atom in a uniform magnetic field 

The diamagnetic hydrogen atom was often used as an example for a quantum 
system whose classical dynamics is chaotic (see e.g. [8] or [9] for an overview). 
As a real physical system it was the topic of studies in experimental physics 
[10,11]. It has even been used for the numerical test of the H correction terms Ci 
and Ci^^ [5,6]. Because of its simple scaling property, which is also fulfilled 
for the h corrections, it is possible to compare the semiclassical results for 
individual orbits with exact quantum mechanical calculations. 

The hydrogen atom in a uniform magnetic field has a continuous symmetry, 
namely the rotational invariance around the magnetic field axis. This symme- 
try can be used to formulate the dynamics in a two-dimensional coordinate 
system. This method works very well for the leading order but it leads to 
new difficulties for the first-order h correction. This aspect is discussed in sec- 
tion 3.1 before wc will look at the "two-dimensional hydrogen atom" with its 
discrete symmetries. 



3.1 The rotational invariance and correction terms 



Written in the cylindrical coordinates g = y/x'^ + y'^, ip and z and atomic 
units, the classical Hamiltonian of the hydrogen atom in a uniform magnetic 
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field in z-direction has the form: 



pI , pI , JL_1+^ 

2^2 r 8 



+ f + 7^-r + ^7V, (47) 

where 



r = y^^^^ (48) 
7 = ^, Sof«2.35-105T. (49) 

-Do 

In this expression is the z component of the angular momentum, which is 
a constant of motion. The paramagnetic term is not considered in equation 
(47) because it is constant. If one exploits the classical scaling property 

~^^2/3- ~*^^-i/3-^ t = 7i, (50) 
the scaled Hamiltonian 

' 2 2 2g^ f 8^ ^ ' 

only depends on the scaled energy e = ^"'^/^E, but not on the energy E or 
the magnetic field strength 7 separately. The regularization of the Coulomb 
singularity can be achieved by introducing semiparabolic coordinates [8] 

ji — Vr + V — \Jf — z (52) 

and the scaled time r 

df = 2f dr = (//^ + i/^) dr. (53) 
This transformation leads to the new Hamiltonian 



(64) 



where the momenta are defined by 



P. = ^, P. = ^. (55) 

The scaled time r is the parameter for the integration of the differential equa- 
tions of the classical values. In the case 1^ — 0, the and v coordinates, which 
are only defined as positive coordinates in equation (52), can be extended to 
negative values. Then, and "py have the same structure as the momenta of 
Cartesian coordinates. This structure of the momenta is very important be- 
cause it is often used in the derivation of the fi corrections (see e.g. equation 
(15) for the monodromy matrix). However, it cannot be found in the quantum 
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mechanical analogue of equation (54). Starting again with the Hamiltonian in 
cylindrical coordinates 

and using the transformations (50) and (52), one arrives at the Schrodinger 
equation for the (//, i/)-part of the wave function: 

f 7^/^ (id d Id d rn^ . , 2\ 

+ ^/iV(//2 + i.2)L(^^^) ^2V'(/i,i^). (57) 



2 djji^ djji ' u du du iJ? z/^ 

1 



In this quantum mechanical expression, the momentum operators for /i and 
u differ from those in Cartesian coordinates even for m = 0, which is the 
quantum analogue of 1^ = 0. As was mentioned above, this difference does 
not have any effect on the leading order of the semiclassical level density. The 
results of the classical calculations agree very well with values which were ex- 
tracted from quantum calculations. But there is no agreement at all between 
the classically evaluated correction terms Ci and Cf^^ and quantum mechan- 
ical calculations. This behaviour can be understood, when one introduces the 
new wave function 

^M = ^, (58) 

which leads to the following Schrodinger equation: 




+ ^[i^v'^{[i^ + i/^)+|u(/x, v) = 2w(/x, v). (59) 

Now, the momentum operators agree with the structure in the classical case 
but an additional term — |7^/^(l//i^ + 1/z^^) appears, which leads to a non- 
vanishing centrifugal part of the potential even for m = 0. The additional 
term has also to be inserted in the classical Hamiltonian: 

(60) 

Since the new centrifugal part is proportional to h"^, which was set to one when 
atomic units were introduced, it contributes to the first-order h correction in 

the expansion of the level density as one can easily see when one repeats 
the calculations from section 2. That is the reason why the influence of the 
continuous symmetry does not appear in the leading order but in the first- 
order corrections. 
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The centrifugal term consists of an integral over the classical periodic or- 
bit with the kernel + l/z^^- which diverges at the coordinate axes. This 
problem does not arise if one regards the hydrogen atom as a purely two- 
dimensional system because in this case there is no centrifugal term. The 
momenta have the structure of the momenta of Cartesian coordinates in the 
classical dynamics as well as in the quantum case. Therefore, in the following 
sections we only look at this two-dimensional hydrogen atom for the calcula- 
tion of the first-order h corrections, as it was done in [5] before. 



3.2 Classical dynamics and Schrddinger equation of the two-dimensional hy- 
drogen atom 



The two-dimensional diamagnetic hydrogen atom can be formulated similar 
to the three-dimensional one. For a magnetic field B = Bex, the classical 
Hamiltonian is given by: 



^ = f + f-^^'V. (61) 

If one uses the classical scaling property (50) and semiparabolic coordinates 
for the two-dimensional case [12] 

= ^g + x, V = ^Jq-x, (62) 
the transformed Hamiltonian reads: 

^=\pI+ \vl - ^/^^ + ^0 + + = 2- (63) 



The Schrodinger equation associated with the classical Hamiltonian (63) can 
be obtained by the same procedure as in the three-dimensional case. The result 
is: 

In this case, the momentum operators have exactly the same structure as for 
Cartesian coordinates. Note that 7^/^ takes the place of fi, which is equal to 
one in atomic units and which is often called "effective W . 
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Fig. 1. Equipotential contours of the potential (65) for different scaled energies e 
are plotted in figure (a). The shadowed area in the coordinate system of figure (b) 
marks the fundamental domain. 

3.3 The symmetry of the potential 

The "potential" of the two-dimensional diamagnetic hydrogen atom in semi- 
parabolic coordinates (see equation (63)) 

V{n, u) = -e {n' + u') + ln'u' {n' + u') (65) 

has a C4t,-symmetry. This symmetry can be seen in figure 1 (a), in which a few 
equipotential contours of the potential (65) are plotted. For the leading order 
(Gutzwiller's trace formula), it is known, how the symmetry of the system can 
be exploited in the calculation of the classical values. In section 4 we look at 
the influence of the discrete Ciu-symmetry on classical values which contribute 
to the first-order h correction terms Ci and Cf 

One has to integrate differential equations along the periodic orbits of the 
classical system in order to calculate the required classical values. It is known 
that, because of the symmetry, the classical calculations can be reduced to a 
fundamental domain [13], which is shown in figure 1 (b). It consists of one 
eighth of the full coordinate plane. If during the integration along an orbit 
one arrives at one of the borders of the fundamental domain, one reflects the 
orbit at the border. This is possible because the borders of the fundamental 
domain are reflection planes of the potential. 

In addition, the symmetry of the potential allows of a compact labelling of 
the unstable classical periodic orbits. In the symmetry reduced fundamental 
domain, the orbits are described by a ternary code, which was introduced by 
Eckhardt and Wintgen [14]. It uses the symbols "0", "+" and "— " and has its 
origin in the description of the orbits of the four disk scattering system, which 
has the same symmetry as the diamagnetic hydrogen atom in semiparabolic 
coordinates. 
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The restriction to the fundamental domain and the introduction of symmetry 
operations lead to new periodic orbits, which are only periodic if one exploits 
the symmetry properties of the system. The periodic continuation can be 
achieved by reflections at the borders of the fundamental domain. An exam- 
ple is the orbit +, which is shown in figure 2 (a). Starting for example on the 
yU-axis the orbit has to be reflected at the angle bisector. After returning to 
the ^-axis a second reflection (this time at the /i-axis) is necessary such that 
the momenta at the initial and final points agree. For practical purposes it 
is often easier not to restrict the calculation to the fundamental domain but 
to find a periodic continuation of the orbit by mapping the final point of the 
orbit on its initial point via a symmetry transformation from C^^, namely ro- 
tations by multiples of 90 degrees (04, C4 = C2, cf), reflections at the coordinate 
axes (ay), and reflections at the angle bisectors {(Jd). Then one can obtain the 
same new periodic orbits as described above. For example figure 2 (b) shows 
the orbit from figure 2 (a) in the case where it is not restricted to the funda- 
mental domain. The periodic continuation is done by a clockwise rotation by 
90 degrees. This method requires only one symmetry operation to render an 
orbit periodic, and is in general easier to implement as the restriction to the 
fundamental domain. 

Figure 2 (c) shows the orbit which corresponds to the double repetition of the 
orbit + in the fundamental domain and figure 2 (d) shows the orbit which 
consists of four times the same repetition in the fundamental domain and 
which is periodic in the plane of semiparabolic coordinates. 



3.4 Symmetry properties of the wave functions and calculation of the quan- 
tum spectra 

Because of the C4^-symmetry of the Hamiltonian (64), the eigenfunctions of 
the system split up into subspaces belonging to a representation of the sym- 
metry group. The symmetry group 6*4^ has four one-dimensional representa- 
tions, namely Ai, A2, Bi and B2, and a two-dimensional representation, which 
is called E. Applying a symmetry element from C^y to a wave function with 
symmetry E leads in general to a linear combination of two (energetically 
degenerate) wave functions. 

At fixed scaled energy e, equation (64) can be considered as a generalized 
eigenvalue problem in the variable 7^/^. The eigenvalues can be calculated by 
diagonalizing a matrix representation of the Hamiltonian in a complete basis 
set with the Lanczos algorithm [15,5]. Using only the wave functions of one of 
the subspaces, which means that the block diagonal form of the Hamiltonian 
is exploited, reduces the dimension of the eigenvalue problem and leads to 
separate spectra for each of the subspaces, which is necessary for the analysis 
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Fig. 2. The periodic orbit + in semiparabolic coordinates for a scaled energy of 
e = —0.1. The dot marks the nucleus at the origin. The orbits shown correspond to 
the single (see figures (a) and (b)), the double (c) and the quadruple (d) repetition 
in the fundamental domain. 

of the influence of the symmetry. 

Up to a maximum 7^^/-^ (the eigenvalues are needed in this form) of about 156 
one obtains circa 12,000 eigenvalues in each of the one-dimensional subspaces 
Ai, A-2, Bi, and B2. In the subspace which belongs to the two-dimensional 
representation E one finds roughly 19,000 eigenvalues up to a maximum 7"^/^ 
of about 140. 



3.5 Analysis of the quantum spectrum by harmonic inversion 



If one looks at the semiclassical level density 



^osc(^) = 



Tr 



Op 



=_{cos(l5f)(To)-|;., 



r 7rn^|det(m^(ro) - 

- hCnesin (^^^(ro) - ^^le) + O (/i^)!, (66) 
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which one obtains from the semiclassical approximation of the quantum Green's 
function (22) via the relation 



-Im (IV G(E)) = g(E) = ^ - ^n), (67) 



one recognizes that it should be possible to extract the amplitude of the leading 
order (Gutzwiller's trace formula), 

T 

^^=|det(m,(ro)-l)|' ^''^ 

(cl) 

and the reduced action S), from a Fourier transformation of the quantum 
level density. A higher precision in both the amplitude and the action can be 
achieved with the harmonic inversion method (see e.g. [9]), where the solution 
of the nonlinear set of equations 

^osc(^n) = E^^e-"^' (69) 



yields the parameters Ag and Si. The Wn are chosen on an equidistant grid. 
For a spectrum which consists of (5-functions such as that in equation (67), one 
has to modify it in order to obtain non-vanishing contributions on the grid, 
e.g., by a convolution of the 5-functions with a Gaussian or by applying a filter 
to the signal. We used the latter method, which is described in [16]. Using the 
harmonic inversion method makes it possible to gain the classical values from 
an analysis of the quantum eigenvalues and to compare these results with 
those from classical calculations. 

The harmonic inversion method can also be applied to extract higher-order h 
corrections 

Cne^Cu + Clr"^ (70) 

for individual orbits in the following way [5]. The quantum level density for- 
mulated in semiparabolic coordinates depends on the variable w = 7""^/^: 



W-2) = E^T— ^ — r-ri — , . .\ '^(^-^n), (7i) 



n 2 



where Ti is the quantum Hamiltonian which belongs to the classical Hamil- 
tonian (63). Using the classical scaling property (63) of the system, the semi- 
classical approximation of the quantum Green's function and the relation (67) 



20 



lead to the level density 



9osc{E) = E - ; <! COS (wSf\To) - 

- ^C,,sin («;5r(To) - ^/z,) + O [(i)' ) }. (72) 



TT 



Combining the formulas for the quantum level density (71) and for its semi- 
classical approximation (72) shows that the harmonic inversion of 



TT 



1 



2 |(^„(;,,i.)|-^(|l + |L)|^„(;,,i.) 
leads to an expansion of the form (69), where 



6{w-Wn) (73) 



Ae = -A^exp 



ITT 



S, = sf\To) 



(74) 



with Ae from equation (68). This method can provide access to the first-order 
h correction, if one subtracts the leading order and multiplies the result by w. 
The final expression is 



TTWr, 



S{W - Wn) 



2 |(^™(/^,'^)|-2^(|^ + |i)|^n(/i,l 

^===— ^ cos (To) - f , 

e ^|det(m,(To)-l)| V 2 J 



(75) 



and the harmonic inversion leads to the amplitude 



j^A exp 



ITT 



Cm, 



(76) 



from which one can easily extract the first-order h correction Cm because all 
other values are known. 



4 Influence of discrete symmetries on the correction terms 



In this section we investigate the influence of the discrete C4^,-symmetry on 
the correction terms Ci and Cf^^. The h corrections are computed for some 
orbits which are not periodic without a symmetry transformation from C4„. 
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4-1 Symmetry transformations in the calculations of the correction terms 



The classical quantities have to be calculated for periodic orbits and, as we 
have seen in section 3, some orbits are only periodic after the application 
of a symmetry operation during or at the end of the integration along the 
classical orbit. It is well known which symmetry operations (reflections or 
rotations of vectors) have to be implemented for the calculation of the phase 
space coordinates but, as was mentioned in section 2, one has to solve a large 
number of additional differential equations for a "new" set of coordinates if 
one wants to calculate the first-order h correction. For example the linearized 
equation of motion (14) has to be solved in order to obtain the monodromy 
matrix. In addition, the derivatives of the coordinates of the orbit and of the 
monodromy matrix with respect to the period are required. Thus we face 
the question of how the symmetry has to be implemented into the correction 
terms. What are the correct transformations for the additional coordinates? 
What is the transformation of the boundary conditions of the classical Green's 
function? 



4.1.1 Correction term Ci 



The solution of the linearized equation of motion (14) for the monodromy 
matrix is essential for the calculation of the classical Green's function by using 
the formulation (13) and by solving equation (19). The monodromy matrix 
is one of the classical values which contribute to Gutzwiller's trace formula, 
and its symmetry behaviour has often been used. However, we look at the 
transformations of its elements because it is an example for all other variables 
which are necessary for the first-order h corrections. All further elements follow 
the same scheme. 



The monodromy matrix can be obtained by solving equation (14) for four 
linear independent column vectors 



Xm — 



( n ^ 

P/iM 
\PuM J 



with initial values which match the condition M{0) = 1. In our system, the 
monodromy matrix has the dimensions 4x4 and the linearized equations of 
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motion in semiparabolic coordinates are: 



— PfiM, 
i'M = PuM, 

PuM = (^e - ^A*^i^^ - i'^^^ fJ^M - + Aii^^) I'M, ^^^^ 

(3 1 \ 

As an example, we look at the anticlockwise rotation by an angle of 90 degrees. 
In this case the symmetry transformation for the coordinates and momenta of 
the orbit leads to 



i'M — PuM, 



(3 1 \ 

2e - 2^^A*^ ~ I^M + (v^Ji + i^A*^) vm, 

/ 3 1 \ 

PuM = (2e - t;^'^^'^ - j^j ^M + (y^iJ^ + i^/i^) IJ'M, 



(78) 



where Jx and P are the new variables in the rotated system. Because of the 
symmetry invariance of the Hamiltonian, one has to find a transformation in 
such a way that 





— PuM, 




t'M 


— PuM, 




PfMM 




3 -2 -2 

2^" 


PuM 




^ -2-2 

2^" 
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In our example this condition leads to the results: 

pM — —T^M, T^M — fJ'M, PuM — —PuM, PuM — P^M- (80) 

The condition is always fulfilled if one uses for the elements of Xm the same 
transformations as for the corresponding phase space coordinates of the orbit. 

Since the matrices A_ and can directly be determined from equation (20), 
there is no need for further symmetry operations in order to obtain the right 
values for these two matrices if the monodromy matrix was calculated with 

^ — * — * 

all symmetry transformations. The same is true for the matrices and B^, 
which follow from equation (17), and for the classical Green's function, which 
follows from (13). This is evident if one recalls the calculations from section 
2.2. Furthermore, if one applies the correct symmetry transformation to the 
elements of the monodromy matrix at every time when one transforms the 
coordinates of the orbit during the integration, all boundary conditions of the 
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classical Green's function are fulfilled because the monodromy matrix is the 
only part that follows from the solution of a differential equation. All other 
values follow from formulas which include already the boundary conditions 
and are not affected by the symmetry transformations during the integration. 

4- 1-2 Correction term Cf^^ 

In addition to the coordinates considered until now, one has to obtain the 
values of the derivatives X'^^^ and iW^"-* for the second correction term 

The equations of motion of the components of X^^^ 

X^'\t, To) = E,jHjk{X{t, To))xi'\t, To) (81) 

have exactly the same structure as the differential equations for the mon- 
odromy matrix, i.e., we already know the symmetry operations which have to 
be implemented here. 

Looking at the differential equations of X^"^^ and X^^\ namely 

xP{t,To) = J:ijHjki{X{t,To))xl'\t,To)xi'\t,To) 

+ j:ijHjk(X(t,To))xj^\t,To), 
xi'\t,To) = E,,i7,,fc,^(X(t,To))x(^)(t,To)x(^'(t,To)XW(t,To) (82) 

+ 3E,,i7,,H(X(t, To))xi'^ {t, To)xf ) (t, Tq) 

+ E,^H,k{X{t,To))xi'\t,To), 

one recognizes that the terms which contain the xj^^ and the xj^^ , respectively, 
have the same form as the equation of motion of the monodromy matrix. 
Therefore, we have to use the same transformation. 

The differential equations of the elements of M^^^ and M^'^\ 

M^j\t,T) = E,fe(i/,,,^(x(t,r))M,(i,r)xW(i,r) + 

Mj^\t,T) = E,k(^H^M^n{X{t,T))Mi,{t,T)Xj^\t,T)Xi'\t,T) 

+ 2H,,UX{t,T))Mlf{t,T)x!^\t,T) + HM^{X{t,T)) 
X Mi^{t,T)X^^\t,T) + H,ui{X{t,T))M^\t,T)) 

(83) 

are derived in the same way as those for, e.g., the X^ by differentiating Hamil- 
ton's equations of motion, and have the same structure. One can easily see 
that the elements of the derivatives iVf*-^-* and M'^'^^ obey the same symmetry 
transformation as the elements of the monodromy matrix. 
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We can summarize the results of this section as follows: Under symmetry op- 
erations all variables which contribute to the first-order h correction transform 
in the same way as the corresponding phase space coordinates of the orbit. 
Whenever a symmetry transformation of the phase space coordinates of the 
orbit is necessary during the integration of a set of differential equations in 
order to obtain a periodic orbit, one has to transform the variables which were 
discussed in this section simultaneously. 



4-2 An example for a symmetry reduced orbit 



As an example for an orbit which is periodic only when a symmetry trans- 
formation is applied, we look at the orbit -|-, which has been introduced in 
section 3.3 and is shown in figures 2 (a) and (b) in the fundamental domain 
and in the full (//,i/)-plane, respectively. 



As a first test for the correct implementation of the symmetry in the correc- 
tion term Ci, the classical Green's function Q{t,t') is considered. In order to 
illustrate the influence of the symmetry transformation, the elements of the 
classical Green's function of the orbit in the fundamental domain (figure 2 
(a)) is shown in figure 3 (a). The Green's function, which in our case is a 2 x 2 
matrix, was rotated in such a way that the Qij (r, r') components represent the 
direction along the classical orbit at the initial point, which is marked by a 
cross in figure 3 (b) . One recognizes the discontinuities at the positions of the 
reflection which appear due to the change in the meaning of the components 
according to the change of the variables by using the symmetry transforma- 
tion. All boundary conditions, which were mentioned in the sections 2.1 and 
2.2, are fulfilled. At the final point of the orbit, all components have the same 
values as at the initial point. The condition Qij{0, r') = ^ij(To, r'), which is ex- 
pected from equation (10), is also fulfilled. The discontinuities in the elements 
Q2j, which appear due to equation (16), are clearly visible. 



Figure 3 (c) shows the classical Green's function for the same orbit, which is 
not restricted to the fundamental domain but closed by a rotation of its final 
point (see figure 3 (d)). Therefore, the only discontinuity appears at the end 
of the orbit. To show this, the first and the last point (after the symmetry 
transformation) of each component are marked with the same symbol. Again, 
the boundary conditions are fulfilled. 
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Fig. 3. The classical Green's function of the orbit + is shown in two cases for 
r' = 0.55 • Tq. The upper two pictures (a) and (b) represent the orbit in the funda- 
mental domain, whereas, the two lower ones (c) and (d) stand for the same orbit 
in the full plane of semiparabolic coordinates, which becomes periodic by the ap- 
plication of a rotation of the final point. In both cases the elements of the Green's 
function are chosen in such a way that the components Qij correspond to the direc- 
tion along the orbit at the initial point, which is marked by a cross in figures (b) 
and (d). In the fundamental domain (figure (a)), where two reflections are required, 
one can clearly see the discontinuities in the elements of the Green's function at the 
two positions at which the orbit is reflected. The discontinuities correspond to the 
change of the meaning of the components. In figure (b) only the last point is affected 
by the symmetry operation (clockwise rotation by 90 degrees). This effect has been 
visualized by marking the last point of each component with the same symbol as 
the first point of the same component. The last point was plotted after the rotation 
has been carried out. In both cases the boundary conditions are fulfilled. 
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Table 1 

Character table of the symmetry group C^y 



L CiaditdtiiJii 


£3 
C- 




f*^ — nr\ 
i'4 — ^2 






Ai 


1 


1 


1 


1 


1 




1 


1 


1 


-1 


-1 


Bi 


1 


-1 


1 


1 


-1 




1 


-1 


1 


-1 


1 


E 


2 





-2 









Jf..3 Influence 0] the symmetry behaviour oj the quantum wave functions and 
comparison of semiclassical and quantum results for the correction terms 



In the different subspaces, which belong to a representation of the symmetry 
group, the individual amplitudes of the orbits can be found with a prefactor, 
which is given by the character x (see table 1) of the element from C^^ which 
provides the periodic continuation of the orbit in the representation of the 
corresponding subspace. The cause of the prefactor is the projection [13] of 
the level density on the subspace belonging to the representation a of the 
group G (in our case G = C^y) with the projection operator 

Pa-^i:Xa{X)D^{X). (84) 
1*^1 xeG 

Here, \G\ stands for the order of the group (|C4y| = 8) and da for the dimension 
of the representation D^- The element Da{X) operates on the final point of the 
orbit with the matrix representation AI(X^^). Only if the symmetry element 
leads to a periodic continuation of the orbit, it will contribute to the 
level density. For the weight of the orbit, one has to take into account that 
the fundamental domain appears eight times in the complete coordinate plane 
and that always eight "identical" orbits appear in the full (/x, z/)-plane. 

As a consequence of the symmetry (eq. (84)), the orbit +, which requires the 
rotation C4 in order to become periodic, should appear in the subspaces Ai 
and A2 with the prefactor pa = I, in Bi and B2 with pa = —1, and in the 
subspace belonging to the representation E it should not appear at all. Table 2 
shows the results of the harmonic inversion of the level density which consists 
only of the eigenvalues of one of the one-dimensional subspaces. The quantum 
values arc compared with classical calculations. The action is independent 
of the subspace, and its classical value is S^^^^ = 0.67746283. The quantum 
amphtudes A^^^^ include information from the Maslov index /^^ (see equation 
(74)). Therefore, they have complex values, and the modulus and phase are 
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Table 2 

The actions sf^^^ and the amphtudes A^^^^ of the orbit + which were extracted 
from the quantum spectrum are compared with their classical counterparts. The 



classical value of the action is S, 



(cl) 



0.67746283, independently of the subspace. 



Because of the projection to the individual subspaces the modulus of the quantum 



mechanical amplitudes corresponds to = = l^^'^'''] and the phase is 



given by arg(^^'^™'*) = — |//^ + aicg{pa) = arg(^^'''''). Since the quantum mechanical 
phases are only determined modulo 2tt, all arguments are shifted into the standard 
intervtil ( — tt.tt]. The + Iims a vMiiisliiiig ami)liludo iu suI)S|")mc(; E. 



subspace 



-,(qm) 



A 



(qm) 



.(cl) 



arg(4''"^)) arg(Af) 



Ai 0.67746323 0.685264 0.68523409 -0.5003-7r -tt/2 

A2 0.67746335 0.685487 0.68523409 -0.5003-7r -tt/2 

Bi 0.67746326 0.685174 0.68523409 0.4997-7r 7r/2 

B2 0.67746325 0.685091 0.68523409 0.4997-7r tt/2 

Table 3 

Amplitude of the orbit ++ in different subspaces. See table 2 for an explanation. 



subspace 



^(qm) 



A 



(qm) 



,(cl) 



arg (^A^^^^^ arg 



Ai 1.35492668 0.64643 0.64647361 0.9991 -tt 

A2 1.35492857 0.64638 0.64647361 0.9985-7r 

Si 1.35492697 0.64652 0.64647361 0.9990-7r 

B2 1.35492700 0.64639 0.64647361 0.9990-7r 

E 1.35492730 2.58530 2.58589444 -0.0012-7r 



TT 
TT 
TT 
TT 





given by 



A 



(qm) 



arg 



\Pa\A 



A 



(cl) 



(qm)\ ^ _^ 

^ 2 



m + arg ipa) = arg (A^ ^) , 



where Pa is the prefactor which appears due to the projection operator P^. In 
table 2 all amplitudes and phases have the expected values. The agreement is 
very good. All arguments are shifted into the standard interval (— 7r,7r]. 

Table 3 shows the same comparison for two repetitions (counted in the fun- 
damental domain) of the orbit +, which are labelled with the symbol +- 1-. 
While the orbit + vanishes in subspace E, the orbit +- 1- can be found in this 
subspace with the expected amplitude prefactor pa — —4:. 

It is important to use the amplitude with the correct prefactor in equation 
(75) for the extraction of the sum Cm of the two first-order correction terms. 
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Table 4 



Correction term Cfn of the orbit + in different subspaccs. The modulus IC^*^™^] and 
phase arg(C^*^™'*) of the analysis of the quantum spectra are compared with the 



classically calculated value Cj 



(cl) 
hi ■ 



subspace 



(cl) 



a 



(qm) 



(qm) 



(cl) 

he 



- 1 



Ai 
B2 



-0.09443001 0.09455 1.0006 • tt 0.0013 

-0.09443001 0.09396 0.9999 • tt -0.0050 

-0.09443001 0.09456 1.0006 • tt 0.0013 

-0.09443001 0.09452 1.0003 ■ tt 0.0010 



If this condition is taken into account, one obtains the correction terms in 
every subspace in which the orbit appears. Table 4 shows the results for the 
orbit + in all one-dimensional subspaces of C4„, which are compared with the 
classically calculated value. 



In tables 5 and 6, the semiclassical results for the correction terms are com- 
pared with values which were extracted from exact quantum calculations. Only 
the eigenvalues of the subspace belonging to the representation Ai were used 
for the harmonic inversion of the quantum spectrum. In this subspace all orbits 
contribute with a prefactor of 1 independently of the symmetry element which 
is required to find the periodic continuation. Quantum mechanical eigenvalues 
up to a maximum w of about 156 were used for the analysis. 



In table 5, classically calculated correction terms are given for some orbits 
which are only periodic with a symmetry transformation. The orbits — I — 

and h have the same shape but are traversed in opposite directions. 

In all other cases mentioned in the tables, the symbol is independent of the 
direction because the orbit is identical with its time reversed counterpart. 
The classical values are compared with the results from the analysis of the 
quantum spectrum in table 6. The agreement of the amplitudes from classical 
and quantum calculations is again very good. In most cases the differences are 
only of the order 10~^. In all three cases in which the difference is larger, the 
actions of these orbits lie close to those of other orbits. This makes the analysis 
of the quantum spectrum difficult. The phases of the quantum mechanically 
calculated correction terms reproduce the correct signs of the classical values 
but the differences are of the order 10~^. A similar behaviour was found in [5] 
for orbits which are periodic without any symmetry transformation. 
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Table 5 

The first-order h corrections for some orbits which are only periodic with a symmetry 
transformation obtained by classical calculations. C^^ is the sum of the correction 
terms Ck and Cff*^ (see equation (70)). 



Symbol 






'-'he 


+ 


-0.09003695 


-0.00439305 


-0.09443001 


++ 


-0.3916016 


0.0299637 


-0.3616379 


0- 


0.0184174 


0.0309192 


0.0493366 


++ — 


-0.578572 


0.063221 


-0.515351 





0.25147 


0.08191 


0.33338 


++ 


-0.92396 


0.16120 


-0.76277 


0-+ 


-0.444747 


0.063019 


-0.381729 


00+ 


2.56347 


0.27548 


2.83895 


0-+-, — + 


-1.93292 


0.29321 


-1.63971 


+++ — 


-2.97331 


0.25335 


-2.71997 


0-++ 


-0.319617 


0.036473 


-0.283144 



Table 6 



Modulus IC^^™"*! and phase arg(C^'^"''') of first-order h corrections obtained by har- 
monic inversion of the quantum spectrum. In the last column the relative difference 

(cl) 

between the classical C^^ and quantum results for the modulus is given. 



■»(qm) 



Symbol 



(cl) 

hi 



(qm) 

hi 



(qm) 



c 



(cl) 

m 



- 1 



+ 


-0.09443001 


0.09455 


1.0006 


vr 


-0.0013 


++ 


-0.3616379 


0.36155 


0.9963 


TT 


-0.0002 


0- 


0.0493366 


0.04974 


0.0293 


TT 


0.0082 


++ — 


-0.515351 


0.51564 


0.9883 


TT 


0.0006 





0.33338 


0.31303 


0.0168 


TT 


-0.0610 


++ 


-0.76277 


0.75772 


0.9731 


TT 


-0.0066 


0-+ 


-0.381729 


0.38192 


0.9861 


TT 


0.0005 


00+ 


2.83895 


2.9894 


0.0529 


TT 


0.0530 


0-+-, 0— 


f -1.63971 


1.6622 


0.8646 


vr 


0.0137 


+++ — 


-2.71997 


2.7453 


0.9726 


TT 


0.0093 


0-++ 


-0.283144 


0.28428 


0.9867 


TT 


0.0040 
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5 Conclusion and Outlook 



In this paper we extended the theory presented in [5] to systems with discrete 
symmetries and the influence of these symmetries on the correction terms was 
discussed in detail. The symmetry transformations presented in this paper 
made possible the calculation of the correction terms for a number of or- 
bits which could not be included without symmetry operations. Nevertheless, 
these orbits, which are not periodic in the plane of semiparabolic coordinates 
without a symmetry operation, contribute to some of the subspaces of the 
quantum spectrum. The classical rcsTilts for the correction terms could be 
compared with values which were extracted from exact quantum calculations. 
An excellent agreement between the results of both methods was found. 

In spite of this success, it must be noted that before semiclassical spectra in- 
cluding first-order h corrections can be calculated over the complete spectrum 
a number of problems still remain to be solved. 

On the one hand, as was already mentioned in [5], and was pointed out in 
section 2, the correction term Ci cannot be calculated in the form presented for 
orbits which have a turning point, but the inclusion of these orbits is essential 
for the h correction of the level density. The correction term Ci diverges for 
orbits with turning points. However, the extraction of the h correction terms 
from quantum spectra with the method presented in section 3.4 leads to results 
which are on the same order of magnitude as the values of orbits without 
turning points. Thus we can assume that the reason lies in an insufficiency of 
the theory, and is not a physical property. 

On the other hand, besides the successfully implemented discrete symmetries, 
physical systems often have a continuous symmetry. For example, it is neces- 
sary to take into account the rotational invariance around the magnetic field 
axis for the complete calculation of all first-order h corrections to the semi- 
classical level density of the three-dimensional diamagnetic hydrogen atom. 
As was mentioned in section 3.1, the problem can be considered as an addi- 
tional centrifugal term in the potential. This term leads to diverging integrals 
if the course of the orbit is not changed. A regularization of these integrals 
suggested in [6] leads to good results for a few individual orbits, however, a 
mathematical justification is lacking. 

Furthermore, for the hydrogen atom it would be interesting to look at the 
w level density, g{w) — — Wn), where w is the scaling parameter 

introduced in section 3.5, in contrast to the "energy" level density (73). In this 
case, there is no need to calculate the quantum mechanical matrix elements 
which in our case are required in equations (73) and (75). This was already 
done before for the leading order (see e.g. [9]). This means that one has to find 
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a scmiclassical expression for a level density with weighting factors, gA{E) — 

En {i^n 



\ [pfi + Pt) "^n/ S{E—En). Eckhardt et al. [17] presented an extension 
of Gutzwiller's theory which allows the calculation of the leading order of the 
level density gA^E). It would be desirable to include higher orders of the h 
expansion in that theory. 
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